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Process of formation of the universe with its further expansion in the first evolution stage is 
investigated in the framework of Friedmann-Robertson- Walker metrics on the basis of quantum 
model, where a new type of matter is introduced, which energy density is dependent on velocity 
of the expansion. It is shown that such an improvement of the model forms potential barrier for 
the flat universe at k — (in contrast with generalized Chaplygin gas model) . Peculiarities of wave 
function are analyzed in details, which is calculated by fully quantum (non-semiclassic) approach, for 
the different barrier regions and stages of evolution. Resonant influence of the initial and boundary 
conditions on the barrier penetrability is shown (in contrast with Vilenkin and Hawking approaches). 
In order to perform a comparative analysis, how much quickly the universe is expanded by different 
models, new quantum definitions of velocity and Hubble function are introduced. These notions 
allow us to study dynamics of evolution of universe in quantum cosmology both in the first stage, 
■ and in later times. 
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I. INTRODUCTION 



Data of astronomical observations suggest on speeding up character of expansion of the present universe. This 
is resent observations of supernova of type la (SNe la) |l|-|3(, more recent cosmic microwave background radiation 
(CMBR) data clusters of galaxies Q, etc. This has been designed as "dark energy" effect In order to 

— — i' explain this phenomenon, different approaches in cosmology have been intensively developed. One can divide them 
into two groups. To the first group we can include models, describing evolution of expansion of the universe at present 
time. However, these models mainly are not quantum, and they do not describe the initial stage in details. One of 
the characteristics used in these models is the parameter of Hubble: from comparison of its calculated values with 
observational data some conclusions are made about the restrictions applied on the considered models (for example, 
see @) and, sometimes, information about character of the evolution of the universe at present times is extracted. To 
C*~) ■ the second group one can include the models oriented on study of the formation of the universe and its further evolution 
in the first stage (i.e., Big Bang). Here, start of expandin g o f the universe is associated with tunneling transition 
through the barrier, usually based on quantum cosmology [lOt [Tlj . In particular, string and brane models 
and also multidimensional models [l8j are intensively investigated in this group. However, in frameworks of such 
models it is usually very difficult to hold out the calculations to the present times, and to connect them with data of 
astronomical observations. So, a bridge between these two research directions is needed. 

The quantum description of the formation of the universe and its evolution in the first stage is a particular difficult 
problem. From the literature we find that the expanded number of papers with practical calculations was performed 
in frameworks of the semiclassical approximation (for example, see [19I423I ]). However, for more accurate study one 
would like to renounce such approximations, and, sooner or later, it will have to do it. Sometimes researchers use 
exactly solvable potentials (for example, see [l9|, I24h29|). But the total number of such initial potentials can be 
counted on the fingers, while proper grounds for such research should be tools for work with barriers of arbitrarily 
shape and without semiclassical approximations. This has caused interest to development of methods of quantum 
mechanics specially oriented on quantum cosmology (for example, see [3(il - l32j |). 

One can understand the most clearly the quantum properties of the formation of the universe in the Friedmann- 
Robertson- Walker (FRW) model. For example, in [33j the resonant behavior of the penetrability of the barrier in 
dependence on the chosen initial and boundary conditions (i.e. coordinate of the start of wave, position of the external 
boundary outside the barrier, energy of the radiation) was opened. This result can strongly (up to several thousand 
of percents) change the results obtained by other authors in the semiclassical approximation. Another result of this 
paper (and [3J,|35|) is that space-time in the first stage of evolution of the universe is discrete rather than continuous. 
But, in later times this property gradually decreases and is lost. At the same time, the semiclassical methods are not 
sensitive to such peculiarity. 

However, as mentioned in [36], in such an approach the barrier is not formed for the flat universe at k = 0, and 
other theoretical grounds should be constructed than quantum mechanics provides this. In this paper, we introduce a 
new type of matter, the energy density of which depends on the velocity of expansion of the universe. We show that 



2 



the inclusion of such a type of matter forms the potential barrier, and such an approach would be an some alternative 
to generalized Chaplygin gas model [13, HH. For a comparative analysis, how fast the universe is expanding in 
frameworks of the different models, we introduce new quantum definitions of the velocity of expansion of the universe 
and the function of Hubble. This basis allows us to investigate dynamics in quantum cosmology, with hope to become 
a basis for testing the quantum models via astronomical observations. 



II. FORMULATION OF QUANTUM DYNAMICS OF EVOLUTION OF UNIVERSE IN 
FRIEDMANN-ROBERTSON- WALKER METRIC 



A. Model with energy density dependent on velocity of expansion 

We shall start from consideration of FRW model. We write action in the form (for example, see [2fJ, Eq.(l), p. 2): 

c f i ( R ^4 p 6a 2 + 6aa + 6fc 

where R is Ricci scalar, a is the scale factor, an overdot denotes a derivative with respect on cosmic time t, k is 
curvature of spatial section which equals to +1, or —1 (for open, flat and closed universe, correspondingly), G is 
Newtonian constant. The energy density in presence of radiation and dark sector can be presented in the form: 

P (a) = Psc s (a) + -gp, PgCg (a) = ( A + ^ j . (2) 

Here, the first item p g Cg (i) represents the generalized Chaplygin gas, including dark energy term A and dark matter 
term B, the second term p ra( j describes energy density of radiation. 

But, from literature we find that practically all calculations of quantum rates of universe evolution in the first stage 
have been performed for the case of k = 1 only (for example, see [1(3, El)- ^ n such a case, potential used in the 
Wheeler-De Witt equation has a barrier, and formation of universe is described as tunneling transition through it. 
At the same time, as it was indicated in [36| , for k — different scenarios of generalized Chaplygin gas model do not 
give the barrier, and this causes questions to applicability of the quantum description of formation of universe and its 
evolution, in general. By such a reason, for such a special case of the flat universe at initial (i.e. for k = 0) we shall 
investigate a new scenario of evolution of the universe where we include a new component of energy density which is 
dependent not only on the scale factor a, but also on its derivative a, i.e. velocity of expansion. 

Hypothesis: We shall introduce into the model the new form of energy (matter) of unknown nature which density is 
dependent on a as 

p vc \(a, a) = po,vei(a) a n , (3) 

where n is some constant parameter. We assume that such a behavior of density can be caused by presence dissipative 
(self-) interacting forces or forces of inversed form (i.e. quantum sources) in this type of matter. In particular, we 
assume that such a type of matter could be some variant of realization of dark energy /matter, and could be connected 
with origin of curvature of the initially flat space-time. In this paper, we shall investigate if inclusion of such a type 
of matter allows us to form barrier for flat universe which produces non-unite penetrability with non-zero reflection. 

Now let us compose the total energy density. Here, we include dark matter and dark energy terms separated 
explicitly, term of radiation and the new density ©: 

/ • \ , Pdust . Pr&d / \ • r> /\, / \ • r> , * 

P (a, a) = pvac + -^j-j + ^j^y + po,vel(a) a = pi (a) + po,vel(aj a . (4) 

Here, p V ac is energy density of vacuum which can be connected with cosmological constant as p mc = A/(8ttG) and 
represents dark energy term, pdust is constant factor of dark matter in dust-like form. Substituting (|3|) to {1}, we 
obtain lagrangian: 

L (a, a) = Af (-a 2 + k- a 2 [p(a) + p ,vci(«) «"] J • (5) 

In particular, for flat universe of closed type (i.e. at p = p V ac, ct = 0, k = 1, J\T = 1, po.vci = 0) we have: 



L (a, a) 
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-a + k — pvac a 
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"-(l-d 2 -Hl c a 2 ), H vac = J± (6) 
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that coincides with eq. (11) in [2(| (up to normalized factor, see p. 4 in that paper). Defining hamiltonian as 
h {a, p a ) = pa - L(a,a) = TV a {-g^j [« 2 + k] + a 2 p(a) + (1 - n) a 2 Po,vel(a) a™ j, 



(7) 



where p is momentum conjugated to generalized coordinate a, we find: 

2ttG 



at n = 1 : h (a,p a ) 
at n = 2 : h (a,p a ) 
3 k 



37V < 



jp 2 + 2 Pa TV a 3 po, ve i(a) + N 2 a 6 pl <vel {a) + k a 



2 37V 2 4 
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47V t 



a 2 /9 ,vci(a)+3/(8^G) 



pi + 47V 2 a 2 a 2 p , V ci(a) + 3/(8tt G) 



(8) 



In particular, hamiltonian of our previous model [33] (i.e. without density p vc i) can be obtained at po,vci — > 0. Here, 
from each expression in (|8]) we obtain: 



/i(a,p a ) 
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Now we apply quantization, after which we obtain stationary Wheeler-De Witt equation (for example, see [20| . 
(16)-(17) p. 4; see (U H^, HTJ] ) . After multiplication on corresponding factors we obtain: 



at n = 1 : 
at n = 2 : 
model [33] : 
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^27Va 3 po,vel(a) ^ +7V 2 a 6 p^ vcl (a) + fc ° 2 ~ ^ a4 /°(a) | <p(a) = 0, 



37V\2 



37V 2 



— + 47V 2 a 2 
oa z 



a 2 Po,vel(a) + 3/(8tt G)] (— - a 2 p(a)) j p(a) = 0, 
127V 2 



(10) 



^ + F ( a) U,(a)=0, F ( a) = - 



8ttG 



fc a + Prad + 



B. Case n = 2 at k = 



Further in this paper, we shall be interesting in the model with density (j4j at k — and n — 2. For the second 
equation in (|10p we have: 



V(o) = -127V 2 



a 2 Po,vei(o) 



+ 



1 



8ttG 



Pdust a + Prad + Pvac 0, 



(11) 



1. Energy density p vo i at large a: vacuum energy vs new type of energy. 

In particular, at large a potential (fTTj) without vacuum energy density p vac obtains the same behavior as potential 
of model [33[ with this vacuum energy term. Comparing potential (lll[) with the third equation in (|10[) . we obtain: 



37V 2 
2ttG 



Pvac a 4 = - 47V 2 a 2 p , ve i(a) p dus t a) ( Po.vd(a) = a) 

' \ 87TG OHuBt / 



' G Pdust 



(12) 



Po,vei(a) — c vo i • a, C vo l 



3 p v 



8ttG pdust 



(13) 



So, at choice of po,vei(a) in form (fT3"]l . the potential (fTTj) without vacuum energy density has the same shape at large 
a as potential of the model [33| . 
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2. Energy density p vc i at small a: formation of the barrier. 

One can form barrier introducing density po,vei in the form: 

dust-like energy: p 0lVe i(a) = p /a 3 or p ve i(a, a) = p a 2 /a 3 , 

radiation-like energy: po,vci(a) = Po/a 4 or Pvci(a, a) = Po a 2 /a 4 , (14) 

other types of matter: po, vc \(a) = p a /a m or p ve \(a, a) = po a 2 /a m , m > 2. 

From here one can write a generalized form of energy density of dust-like matter or radiation, which forms barrier at 
small a: 

Pdu™t( a > ") = Pdust [l + Cd us t • a 2 ] /a 3 or pf^ (a, a) = p rad [l + c rad • a 2 ] /a 4 . (15) 

Now we compose density po vei; adding its behavior (|13[) at large a: 

Po,voi(a) = Pdust Cdust/a 3 + Cvci • a or po, ve l(a) = Prad c rad /a 4 + c vc i ■ a- (16) 

In particular, using the first formula for a generalized dust-like matter, we write potential (fTTj) (at J\T = 1, 8ttG = 1, 
Pvac = 0): 



7(a) = -12 |p rad + Pdust a| 1 + p 0: vci(a) • a 2 /3 = -12 j 



Prad Pdust Cdust 
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i Pdust c dust 
Prad H ^ 



, Prad C V cl 3 , Pdust Cvci 

Pdust a H a H 



a'} 



3. Comparison with Chaplygin gas model and choice of parameters 

Some estimations of rates were performed on the basis of analysis of dynamic of wave packet tunneling through the 
potential barrier [TTj] . However, a case of a = 1 (i.e. for standard formulation of Chaplygin gas) corresponds to those 
calculations (see energy density (14) in that paper). So, we shall use parameters of that paper for our calculations, 
in order to give other researchers more material for comparison and analysis. According to that paper, let us rewrite 
energy density of the generalized Chaplygin gas in the following form: 

/ d \ 1/(1+") i / R \ 1/(1+") 

**<■> = (* + 371*0 =^ + ^) • < 18 » 

where A = Air 1+a and B = Bir 1+a . We have: 

iw 11 ( \ , Bl,(1+a) Pdust 

1) at small a: p gCg (a) -> = —3-, ^ ig ^ 

2) at large a: p gCg -> ee p vac . 

In [ll| the following values A = 0.001, B = 0.001 were chosen (see Fig. 1 in that paper). From here we find 
parameters: 



(20) 



p dust = s 1 /(i+«) = s 1 /(i+«) 7r -(i+«)-i/(i+«) = fli/(i+«)/ ff = o.001 1 /(i+«)/ T , 

p vac = AV(l+«) = A 1 /(l+«) jr-d+aJ-l/Cl+a) = ^l/(l+a) /?r = O.OOI^^+^/tt. 

C. Operators of the function of Hubble and velocity of expansion of the universe 

Let us consider how the velocity of the evolution of the universe can be defined in quantum approach. In classical 
mechanics the velocity of the particle is related to its momentum. In quantum mechanics, there is connection between 
the corresponding operators. According to basic positions of quantum mechanics, determination of the wave function 
\& at some moment of time to, not only fully describes all quantum properties of the studied system at this time to, 
but also fully determines its evolution in all future times. In other words, derivative /dt of the wave function on 
time at any given moment of time to is determined by this wave function ^ at to- If w e want to use the principle of 
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superposition, this dependence should be linear. On this basis, instead of the stationary Wheeler-De Witt equation 
we shall be based on the following equation: 

inf^L = h% (21) 

at 

where h is some unknown operator. Let us clarify what this operator should be to. If to assume that the universe 
expands classically (i.e. with high degree of confidence we can neglect by the quantum properties) at large values of 
the scale factor a, then the wave function can be written as 

»f = Ae lS/h , (22) 

where S is action. Substituting this expression for the wave function into (|21[) and neglecting by change of amplitude 
A over time, we find: 

d JL = i d A % (23 ) 

dt h dt ' y ' 

Comparing this expression with definition (121[) . we conclude that in the limiting case the operator h is reduced to 
simple multiplication on value of —dS/dt, i.e. Hamiltonian function. Now we have introduced time into quantum 
equation, connecting it with Hamiltonian operator and taking into account that the asymptotic representation of the 
wave function is connected with action as (|2"2")l . 

As the non-stationary quantum equation has already been defined, now we can define operator of the velocity using 
general rule of differentiation of operators over time (for example, see (19.1) in [4l|, p. 78-79). 



Definition 1. We define operator of the velocity as 



jr(ha-ah). (24) 



Here, we use requirement (as in standard quantum mechanics): 

^- t J\^{a,t)\ 2 da = Q. (25) 



Substituting hamiltonian in form (|5J) or ([§]), we find 

i 8ttG f d 



at n= 1 : a = ^ ^-f- i + iAfa 3 po,vei(a) 
n (xN a I da 



at n = 2 : a = 

h 2Ma 



i 1 d 



a 2 p ,vei(a)+3/(87rG) 



model [331 



i 8nG d 



a, = — 



U QN a da 

Definition 2. We define operator of the function of Hubble as 



H(a) = -a. (27) 



a 



According to this definition, we shall consider the parameter of Hubble at the given scale factor eto, as action of certain 
operator H on the wave function at ao- The wave function is not eigenfunction of operators of the velocity and the 
function of Hubble, as there are no any constant eigenvalues for these operators. So, action of these operators on the 
wave function can be written as 

a * (a) = v(a) * (a), H (a) * (a) = H (a) * (a). (28) 

Here, v(a) and H (a) are some functions changed in depending on a. Near to arbitrarily chosen value ao these functions 
tend to certain well-defined fixed values, which can be locally considered as eigenvalues of operators of the velocity 
and the function of Hubble at ao. Thus, on the basis of the functions v(a) and H (a) we shall understand the velocity 
and the parameter of Hubble in the quantum approach. For practical calculations, one can obtain these functions as 



(29) 
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D. Quantum definition of duration of existence of universe 

As the function of Hubble has been defined, as next step we can define duration of existence of the universe. We 
shall be interesting in such a characteristic which is dependent on the scale factor a: this should allow to see clearly 
behavior of expansion (evolution) and compare such a dynamic for different model scenarios. So, we introduce the 
following characteristic: 

a 

*v>=Jim«- (30) 

One can see that such a definition is analog of the classical definition of the universe age given in classical cosmology 
(for example, see (36), p. 11 in (42j). 

E. Rescaling 

Inside region < a < 100 potentials of the considered models above achieve essential values (that can cause serious 
difficulties in practical calculations and further analysis). But, they can be decreased via renormalization of the 
corresponding equations (|10[) (we shall call such a procedure as rescaling). By such a reason, let us pass to a new 
variable: 



'12 AT 2 

«new = Vd \d, V 



Now we obtain new equations for determination of the wave function: 

& ( 87T G \ f . n . - / Anew \ ^ / a ncw \ d , r n / fl ncw \ ^ 



(31) 



at n = 1 



d 2 „ rl fa 

at n = 2 : - — h 8vr G 



model [33] 



)j ■ r, \r ( a ncw \ 3 / "new \ 9 , \ r2 f "new \ 6 2 ( "new \ 

f 3A/"\ 2 , / Anew \ 3 A/" 2 / ^new \ / ^new \ 1 / \ r\ 
+ UTgJ *{—) -2^Gl^rJ p(— )j^(a now ) = 0, 

i(^f)%o,ve 1 (^)+l/(87rG)" x 

/ ^new \ ^ / Q'now \ ^ / Q'ticw M / \ r\ 

KIT) "l^rJ H— JW a new)=0, 

^^ 2 _ | ^ (flnow) = Q 





r 3fc 


X 


-8ttG 


d 2 
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8ttG 

For operator of velocity (|26[) . determined via new variable, we obtain: 
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(32) 



at n = 1 : d 
at n — 2 : d 
model [33]: d 



7i 6A/"a ne w 1 9a 
i v 



new ^ L'uaiew 

2 r v 2 



■ » /* / ^ncw \ / ^new \ 1 

+ XHV^ — ) P0.vcl( — J , 



?i 2A"a nG w 
i 8ttGi/ 2 9 



(— A).vci(— ) +3/(87rG)l (33) 



h 6A"a ncw da ncw 

III. THE MODEL [32 

A. Tunneling deeply under the barrier 

Before applying formalism above to the model with energy density dependent on velocity, at first we shall clarify, 
how this formalism is working in the model (33|. In the beginning, let us consider tunneling which take plase 
deeply under the barrier. We shall choose value for energy radiation as £7 ra d = 100. Results of such calculations of 
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FIG. 1: (Color online) Penetrability of the barrier for model [33] at -E ra d = 100 (parameters of calculation: 10000 intervals 
at a max = 100): (a) dependence of penetrability on the starting point a star t in region from zero up to internal turning point, 
at fixed a max = 100: at increasing of a star t the penetrability oscillates, maxima are increased and minima are decreased, (b) 
dependence of the penetrability on the starting point a star t in more detailed consideration, (c) dependence of the penetrability 
on external boundary a max at fixed Ostart =0.1: one can see that at increasing of a max the calculated penetrability tends to 
some definite limit value, which we chose for further calculations and analysis. 



penetrability are presented in Fig. [TJ In the first figure (a) it is shown how the penetrability is changed in dependence 
on the starting point a s t ar t- One can see that it has oscillating behavior, maxima are slowly increased and minima 
are decreased at tending of the starting point to internal turning point (at fixed a max = 100). The second figure 
(b) is the same result, but in more enlarged presentation. In the last figure (c) it is shown how the penetrability is 
changed at increasing of the external boundary a max , starting from the external turning point (at fixed a s tart =0.1, 
corresponding to coordinate of minimum of internal well before the barrier). From here it follows that: (1) our 
approach gives convergent calculations, (2) determination of the penetrability on the basis of shape of the barrier 
inside tunneling region only (like semiclassical approaches of the second order) seems to be enough far from calculated 
penetrability, obtained after taking into account external tail and internal shape of the well, (3) inspite of sharp 
decreasing of the potential in the external region (after external turning point) at increasing of the scale factor a 
(which decreases with evident acceleration!) calculations of penetrability are convergent and they coincide to some 
definite value (that allows us to speak about reliable values of the penetrability for such potentials). 

Now we shall analyze the function of Hubble defined by formula (|29p . In Fig. [5] it is shown how the modulus of 
this function is changed in dependence on scale factor a. In the first figure (a) general picture is shown: at small 




Scale factor, a Scale factor, a Scale factor, a 

FIG. 2: (Color online) The modulus of the function of Hubble in dependence on scale factor a at i? ra d = 100 (calculation 
parameters: starting point a star t = 0.1; 10000 intervals at a max = 100): (a) whole range of variable a can be separated on 
3 regions: in the internal region (at a < 21) presence of chaotic peaks and minima is observed, in the tunneling region (at 
21 < a < 55.6) the function is smooth and has no any oscillation, in the external region it increases monotonously (slowly 
transforming to linear dependence) , (b) at transition from the tunneling region to external one smooth minimum of this function 
is observed, corresponding to external turning point, (c) location of minima and peaks at small a is similar (the last minimum 
corresponds to the internal turning point). 

a (close to a = 21) sharp chaotic peaks are observed, then this function has one clear maximum and minimum, at 
finishing it increases monotonously. In the first consideration, such a behavior of the function of Hubble (especially 
at small a) looks to be enough strange... But after increasing insight of the last right minimum at a — 55.6 one 
can see (see figure (b)) that it is stable and is not equal to zero (that confirms coincidence and stability of computer 
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calculations), corresponding to the external turning point — i.e. it separates the tunneling region from the external 
region! In the external region the modulus of the function of Hubble increases monotonously, slowly transforming to 
linear dependence. Inside the region 21 < a < 55.6 the modulus of the function of Hubble has no any oscillation — this 
is the tunneling region, which is finished by convergent minima at both sides. Detailed analysis of the internal region 
(at a < 21) shows that location of minima and peaks here is similar (see figure (c)), the last minimum corresponds to 
the internal turning point. 

Here the following question can be appeared: whether peaks are finite at small a or we are dealing with divergencies 
in calculations? Let us consider formula (|29[) for the function of Hubble: one can see that maximal values should be 
caused by practically zero values of the wave function, which is in the denominator. The wave function calculated 
for this process is shown in Fig. [3J In the first figure (a) the modulus of the wave function is shown: indeed, with 




Scale factor, a Scale factor, a Scale factor, a 



FIG. 3: (Color online) The wave function in dependence on the scale factor a at i5 ra d = 100 (calculation parameters: the 
starting point a s tart = 0.1; 10000 intervals at a max = 100) (a) in the internal region (at a < 21) accelerated growth of the 
maxima of modulus of the wave function (with minima very closed to zero) up to the internal turning point is observed, in 
the tunneling and external regions (at 21 < a) this modulus decreases without oscillations, (b) in the internal region the real 
and imaginary parts of the wave functions oscillate almost simultaneously, while in the tunneling and external regions they 
extremely fall down, (c) consideration of the wave function in the logarithmic scale clearly demonstrates its non-zero values in 
the tunneling and the external regions. 

clearly determined finite maxima (which should definitely give not zero minima of the function of Hubble) the sharp 
minima are seen tending to zero — it explains presence of the sharp peaks in the Hubble function in Fig. [TJ This 
could be explained by almost simultaneous zeroing of the real and imaginary parts of the wave function. However, 
this is strange as we have non-zero complex wave function (as it defines non-zero constant flux directed outside inside 
entire region of variable a), so zeros of its real and imaginary parts should not be coincide anywhere. 

In the next figure (b) the real and imaginary parts of the wave function are shown. Here, one can see clear stable 
curves that demonstrates convergence and stability of calculations, and we need to understand the result. However, 
the curves behave almost similarly: their maxima, minima and zeroes are located at close coordinates. This situation 
is similar to the behavior of the wave function of the bound state for a particle inside a potential well (with infinitely 
high boundaries). Indeed, at the chosen energy the penetrability is very small (Tb ar ~ 10~ 255 ) and output is extremely 
small. By other words, we are dealing with the quasi-stationary state with extremely small output outside, which is 
very close to stationary one, practically. 

In the last figure (c) the modulus of the wave function is shown in the logarithmic scale. Here, one can see that 
it uniformly decreases inside the tunneling region. It is clear that in the external region the modulus is constant 
and non-zero: this confirms the convergence of the calculations, proves presence of the stable non-zero minima of the 
modulus of the wave function and the finite maxima of the peaks of the function of Hubble in the internal region. 

B. Tunneling near the barrier maximum and above-barrier propagation processes 

Now let us consider a case where the tunneling occurs near the barrier maximum. In this case, we can use 
our previous analysis in [33| and choose i? ra( j = 220. So, at start in point a star t = 0.1 we obtain Tb ar = 
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1.52129237224042 cdotl0~ 7 , i? bar = 0.999999847870763 and the condition T bar + i? bar = 1 holds up to 14 digits 1 . 

Our calculations show that the penetrability in dependence on the starting point a sta rt and on the external boundary 
flmax behaves like to the case studied above at B ra( j = 100 (also see fH}). Placing the starting point in the minimum 
of the internal well, we calculate the wave function and the function of Hubble (see Fig. 2]). From these figures it 
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FIG. 4: (Color online) The wave function and the function of Hubble in dependence on the scale factor a at S rac j = 220 
(calculation parameters: the starting point a s tart = 0.1, 10000 intervals at a max = 100): (a) in the internal region number of 
oscillations of the wave function is much larger and its maxima increase stronger in comparison with results at _E ra d = 100 (see 
Fig. |3](b)), that is explained by increasing of the energy i5 ra d and enlarging of the internal well region; (b) in the tunneling and 
external regions the modulus of the wave function is changed essentially weaker in comparison with results at -E ra d = 100 (see 
Fig. [3] (c) ) , that indicates on essential oncoming of energy -E ra d to the barrier maximum; (c) the function of Hubble behaves 
similarly to the case at i? ra d = 100 (see Fig. [2}, but here the tunneling region is much smaller and the function of Hubble in it 
is smaller. 



is clear that these functions behave as obtained above at E ra d = 100. However, in this case we observe: (1) in the 
internal region number of oscillations of the wave function is essentially larger and the difference between its maxima 
is reinforced (see Fig. @] (a) in comparison with Fig. El (b)), caused by enlarging of the internal region with shift of 
the internal turning point to the right; (2) in the tunneling and external regions the difference between the maxima 
and minima of modulus of the wave function is much smaller (see Fig. (b) in comparison with Fig. [3] (c)), that 
indicates on essential oncoming of the energy -E^ad to the barrier maximum and much stronger outgoing flux through 
the barrier outside; (3) the tunneling region is less, which can be easily found by typical behavior of the function of 
Hubble (see Fig. 1(c)). 

Now let us find out how these pictures will be changed, if we increase the energy -E ra d above the barrier. Results 
such calculations are shown in Fig. [SI where we have chosen i? ra d = 250. In the first figure (a) the function of Hubble 
is shown. One can see that (1) In the internal region the sharp peaks and minima have disappeared completely, 
(2) In the middle region the behavior of this function, typical for the tunneling region, has disappeared also, and 
instead there is a smooth minimum in the coordinate of the potential barrier maximum, (3) In the external region the 
monotonic increasing of the function of Hubble remains, slowly becoming to linear dependence. The wave function and 
its modulus are shown in next two figures (b, c). In this case, we see oscillatory behavior, typical for the above-barrier 
energies. However, its maximum is located at coordinate of the barrier maximum, with a monotonic decrease to both 
sides. 

Now let us analyze, how velocity defined by the formula (|2"9")l behaves in each case. Such calculations are presented 
in Fig. |6] In general, the behavior of the velocity looks like the function of Hubble. From the figures it is clear 
that if in the internal region there are sharp peaks and minima at the sub-barrier energies (see figure (a)), then they 
disappear at increasing of the energy i? ra d above the barrier (see figure (b)). However, after detailed analysis one can 
see the oscillations with a general tendency to decreasing (and decreasing amplitude with increasing i? ra d)- They can 
be explained by wave nature, which quantum-mechanical treatment of the process gives us. They are the previously 
considered sharp peaks and minima at sub-barrier energies, i.e. now picture of sub-barrier and above-barrier processes 
becomes unite. Figure (c) shows that the function of Hubble at different energies -E ra d tend to the same limit behavior 
at large a. 



1 Note that the semiclassical methods usually do not give the reflection coefficient i?b ar an d the mentioned test is not applied. However, 
in discussions on comparison between fully quantum calculations and semiclassical ones this important point is usually ignored, with 
assumptions on advantage of the semiclassical apparatus without alternatives. 
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FIG. 5: (Color online) The function of Hubble and wave function in dependence on the scale factor a at -E ra d = 250 (calculation 
parameters: the starting point a star t = 0.1, 10000 intervals at a max = 100): (a) the function of Hubble has a smooth form: in 
the internal region the sharp peaks and minima (typical for above-barrier energies and reflecting bound state) are disappeared 
completely, in the middle region there is only one smooth minimum in the coordinate of the barrier maximum (tunneling 
maximum is disappeared as shown in Figs. [2] (a) and|4](c)), in the external region the accelerated increase of the function of 
Hubble remains (as in Figs. [2] (a) and|4](c)); (b) the wave function is oscillatory in the entire region; (c) modulus of the wave 
function has maximum in coordinate of the barrier maximum. 




Scale factor, a Scale factor, a Scale factor, a 

FIG. 6: (Color online) The velocity of expansion of universe in dependence on the scale factor a (calculation parameters: the 
starting point a star t = 0.1, 10000 intervals at a max = 100): (a) in the internal region at sub-barrier energies there are sharp 
peaks and minima, in the tunneling region the velocity is smooth, in the external region velocities at different energies tend 
to the same limit, (b) for the above-barrier energy i? ra d = 250 sharp peaks and minima are disappeared, hump of tunneling 
(previously observed for the sub-barrier energies) transforms to smooth minimum in coordinate, corresponding to the potential 
barrier maximum, (c) the functions of Hubble at different energies i? ra d tend to the same limit at large a. 



So, there is the following picture at the above-barrier energy E Ta< i — 250. We have no any classically forbidden 
region inside whole area of the scale factor a. Hence it would seem that this situation forbids to consider formation 
of the universe, and we would reject this case. However, we see that for values of a less than coordinate of the barrier 
maximum (we denote it as abar) the velocity gradually decreases to minimum with increasing of a. In this case, the 
modulus of the wave function increases up to maximum: it points to a gradual increase of probability of appearance 
of the universe, with maximum at abar- If to take the spherically symmetric picture of the extension into account 
(where the density of matter in filled volume should not increase and it can be associated with the probability of 
appearance of the the universe), then such an increase of the probability is contrary to natural expansion of the 
universe in classical treatment. So, this situation is more suitable to description of formation of the universe with 
classical space-time, with its birth at abar- Starting from abar, the expansion of space-time becomes classic. Before to 
this coordinate, the universe is formed with a gradual damping expansion, at maximum in point abar this extension 
is practically stops for a certain period, then further expansion begins with the acceleration. This logic points us 
to competence of quantum description of formation of the universe at the above-barrier energies (when there is no 
tunneling) . 

Also another property is found: All velocities at different energies E ra d tend to the same limit with increasing a 
in the external region: i.e. this model gives the same dynamics of the accelerated expansion of the universe for later 
times with completely different scenarios of its evolution in the first stage. 
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IV. MODEL WITH ENERGY DENSITY DEPENDENT ON VELOCITY OF EXPANSION 



Now we shall analyze the behavior of the model with energy density dependent on the velocity at k = and n = 2. 
The potentials generated in this approach are shown in Fig. [7] one can see that the energy density depending on the 
velocity forms the barrier for different values of the potential parameters. In next Figs.|5]the results of our calculations 
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FIG. 7: (Color online) Inclusion of the density component dependent on the velocity into the model allows to form the barrier 
at k — 0: One can see that with increase of parameter Cd us t the barrier is shifted to the right down, forming more rounded 
shape (calculation parameters: -B ra d = 100, solid orange line is for Cdust = 10 3 , dash-dotted blue line for Cd us t = 10 5 dash-double 
dotted brown line for c dus t = 10 6 , dashed green line for c dus t = 10 7 short dashed purple line for c dus t = 5 ■ 10 7 short dash-dotted 
red line for Cdust = 10 s ). 



of the function of Hubble, velocity and the modulus of the wave function are presented. One can see that the universe 




Scale factor, a Scale factor, a Scale factor, a 

FIG. 8: (Color online) The function of Hubble, velocity and the wave function in dependence on the scale factor a at E la .d = 100 
and k = for the model with velocity-dependent density (calculation parameters: the starting point Ostart = 0.01, 10000 
intervals at a max = 100; blue solid line 1 is for Cd us t = 1, brown dashed line 2 for Cd us t = 10 5 , green shot dash-dotted line 3 
for Cdust = 10 7 , red dash-double dotted line 4 for Cd us t = 10 10 ): (a) the function of Hubble decreases monotonically at increase 
of the scale factor a and parameter Cd us t (curve 5 is added on this figure for the model which is much higher and tends 
to higher limit); (b) for different values of Cd us t the velocity tends to the same limit (at Cd us t = 10 10 this coincidence occurs at 
much higher a), (c) the modulus of the wave function is similar to the barrier shape, maximum of this modulus corresponds to 
the coordinate of the barrier maximum. 



is slower expanding in frameworks of such a model, with increase of the scale factor a its velocity decreases, tending 
to some limit monotonic dependence. So, the model with the density dependent on the velocity gives slow expansion 
at late times, in contrast to the model [331 ] - 

In Fig. IHlthe velocity of expansion of the universe for this model in dependence on the scale factor a at variations 
of parameters c vo i, n and -E ra d is shown. In next Figs. [10] we demonstrate comparing calculations for the velocity of 
evolution of the expansion of the universe by this model with the included cosmological constant (at k = 0). Our 
results of calculation of the duration of the universe defined by (1301) with different parameters are presented in last 
Fig. QH 
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FIG. 9: (Color online) The velocity of expansion of the universe in dependence on the scale factor o for the model with component 
of the energy density dependent on the velocity at k = and variations of parameters c ve l, n and E la ,d (calculation parameters: 
the starting point a s t a rt = 0.01, 10000 intervals at a max = 100): (a) calculations for different c vc i at fixed Cd us t = 10 s , E rad = 100 
and n — 1 (blue solid line 1 is for c vc i = 10~ 5 , brown dash-dotted line 2 for c ve i = 1, green dashed line 3 for c vo i = 10 5 ); (b) 
calculations for different powers n at fixed Cd us t = 10 J , c vc i = 3 and -E ra d = 100 (purple dash-double dotted line 1 is for n = 1/2, 
blue solid line 2 for n — 1, orange dash-dotted line 3 for n — 2, red dashed line 4 for n = 3, green short dotted line 5 for 
n = 4, wine short-dashed line 6 for n = 8); (c) calculations for different energies -E ra d = 100 at fixed Cdust = 10 5 , c vc i = 3 and 
n = 1 (red dashed line 1 is for -B ra d = 1000, blue solid line 2 for E la _d = 100, orange dash-dotted line 3 for -E ra d = 10, green 
short-dotted line 4 for -E ra d = 1, wine dash-double dotted line 5 for i? ra d = 0.01). 
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FIG. 10: (Color online) The evolution of the expansion of the universe for the model with the component of energy density 
dependent on the velocity with the included cosmological constant at k = 0, Cd us t = 10~ 5 and c vc i = 10 -5 (calculation 
parameters: the starting point a s t ar t = 0.01, 10000 intervals at a ma x = 100; orange dash-double dotted line 1 is for -E ra d = 10000, 
blue solid line 2 for -E ra d = 1000, red dashed line 3 for _B ra d = 100, green dash-dotted line 4 for -B ra d = 10): (a) velocity: varying 
the energy £" ra d, one can obtain a smooth transition from the accelerated expansion to decelerated one for later times; (b) the 
potential: one can see a presence of the barrier for different parameters. 



V. CONCLUSIONS 



Process of formation of the universe with its further expansion in the first evolution stage is investigated in the 
framework of Friedmann-Robertson- Walker metrics on the basis of quantum model, where a new type of matter is 
introduced, which energy density is dependent on velocity of the expansion. It is shown that such an improvement of 
the model forms potential barrier for the flat universe at k = 0. Peculiarities of wave function are analyzed in details, 
which is calculated by fully quantum (non-semiclassical) approach. Resonant influence of the initial and boundary 
conditions on the barrier penetrability is observed. In order to analyze dynamics of evolution of the universe in 
quantum approach, we introduce new operators of velocity of expansion and the function of Hubble. This basis allows 
us to study dynamics of evolution of universe in quantum cosmology both in the first stage, and in later times. 
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